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We present results for the light quark masses and the neutral kaon mixing parameter Bk 
from lattice QCD. Our data set includes lighter than physical light quark masses and 5 lattice 
spacings so that chiral extrapolation is not necessary and cutoff effects are fully under control. 
We obtain fully nonperturbative predictions for m u d = (m u + md)/2, m s and Bk in the 
RI scheme with M*, Mk and Mq as the only input quantities. Using perturbative 4- loop 
respectively 2-loop running and dispersive input from r\ — > 3ty, we obtain mjf s (2 GeV) = 
2.17(4)(10) MeV, mf^(2 GeV) = 4.79(7)(12) MeV, mf (2 GeV) = 95.5(1. 1)(1.5) MeV and 
Bk = 0.773(8) (8) where the first error is statistical and the second systematic. 

1 Introduction 

Lattice QCD is a tool to perform ab-initio calculations of QCD in the nonperturbative regime. 
One can stochastically perform the functional integral over gauge and fermion fields on a lattice 
regulated theory in finite volume. In order to obtain QCD predictions, it is then necessary to 
(a) remove the cutoff (b) extrapolate to infinite volume (c) tune the bare parameters of the 
theory or interpolate/extrapolate to these bare parameter values such that a predefined set 
of experimentally accessible, dimensionless quantities (e.g. hadron mass rations) assume their 
physical value. If these requirements are met, one can use lattice QCD to obtain ab-initio QCD 
predictions of quantities not used as input in step (c) with a statistical error that arises from 
the stochastic integration. 

Generally speaking, the challenge for lattice QCD is to simultaneously fulfill all the require- 
ments (a)-(c) with control over systematic errors arising from each step and to minimize the total 
(statistical plus systematic) uncertainty on a target quantity with given computer resources. It 
has been shown recently, that relatively straightforward quantities such as the ground state light 
hadron spectrum can be reproduced with a few percent accuracy^ Here we present the results 
of a determination of the light and strange quark masses as well as for the neutral kaon mixing 
parameter Bk with controlled errors on the percent level. For the full technical details, we refer 
the reader to the original publicationsHMIH 

2 Quark masses 

Light quark masses are fundamental parameters of the standard model Lagrangian that are inac- 
cessible by direct experiment. In order to compute them, we compute some light hadron masses 
in lattice QCD with a number of bare input light quark masses. We then interpolate the bare 
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Figure 1: Continuum extrapolation of the average up/down quark mass, of the strange quark mass and of the 

ratio of the two. 

quark masses at one value of the bare coupling g to the point where the measured light hadron 
masses take on their physical value and renormalize them. We compute the renormalization 
constant nonperturbatively in the RI-MOM schemel^H 

For technical reasons we work in the isospin limit and correct for the small isospin breaking 
at a later stage. We therefore take as input for finding the physical point the ratios of isospin 
averaged hadron masses M n /M^ and Mk/Mq. The lattice cutoff or, equivalently, the lattice 
spacing a at the given bare coupling g itself is determined by comparing the dimensionless mass 
of the f2 baryon as measured on the lattice to the physical mass Mq. We interpolate to the 
physical point with various functional formsP The resulting spread enters (as a subdominant 
part) into our systematic error. Repeating this procedure at various different values of the bare 
coupling (in our case we use 5), we get the renormalized quark mass as a function of the lattice 
spacing a that we can extrapolate to a = (see fig. [I]). 

For our lattice action, the leading term in the continuum extrapolation is formally of 0(a s a). 
There is however strong evidence that this term is numerically subdominant (for lattice spacings 
we work at) to the term of 0(a 2 )PWe therefore use both forms in our analysis and include the 
resulting spread into our systematic error. 

Finite volume corrections to stable particle masses are generically exponentially suppressed 
in M W L and can be corrected for systematically^ (see fig. [2]). We included these corrections and 
found them generically to be at the permil level. 

In order to obtain individual masses for the up and the down quark, we use the double ratio 
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to convert our precise result for m s /m u <i = 27.53(20)(8) into an estimate of m^/m^. In principle, 
Q can be determined experimentally from r/ — > 3ir decays via dispersion relations. Due to 
the imperfect experimental data, there is some amount of modeling involved and we use a 
conservative estimate Q = 22.3(8) from a recent reviewP 

As a last step, we convert the individual quark masses from the nonperturbative RI-MOM 
scheme into the MS scheme. In both schemes, the running of the quark mass is known to 4- loop 
orderHE2] As demonstrated in fig. [3J the nonperturbative running is well described by 4-loop 
perturbation theory above \i = 4 GeV. We therefore compute our quark masses in the RI-MOM 
scheme at \x = 4 GeV and further convert these numbers into the MS scheme using the results 
of Chetyrkin and Retey?^ Our final numbers are 
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Figure 3: Nonperturbative over 4-loop perturbative running of the inverse of the mass renormalization constant 
Zs in the RI-MOM scheme. Above 4 GeV agreement is reached within the statistical precision of our data. 
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Figure 4: Contribution of the various operators towards the final value of Bk for the ensemble with the lightest 
pion mass. The contribution from the standard model operator Oi is clearly dominant. 



3 Neutral kaon mixing 

Neutral kaon mixing is responsible for indirect CP violation in K — > 2tt decays. It is phe- 
nomenologically described by the parameter e, which contains the hadronic matrix element of 
the standard model AS = 2 operator AS=2 = (s7^c?)(s7 m l^) that is usually parameterized as 
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A precise determination of Bk together with an experimental measurement of e thus constitutes 
a precision test of the standard model in the kaon system which is particularly relevant for 
constraining various standard model extensions!. 11 ! 12 ! 1 ^ ! 

We have performed a lattice determination of Bk using the same setup as for our quark 
mass determination^ One particular point to note is that our lattice discretized fermion action 
does only exhibit approximate chiral symmetry that gets fully restored in the continuum limit 
only. Consequently, mixing of the standard model operator, which has the structure 0\ = 
(V —A)(V —A), with other dimension-6 operators that is forbidden in the continuum is allowed at 
finite lattice spacing. These other operators are Oi = VV — AA, O3/4 = SS^fPP and O5 = TT. 
As the standard model operator is chirally suppressed, these mixings can in principle be very 
large. Due to the good approximate chiral symmetry of our actional the mixing contributions 
to Bk are actually tiny as displayed in fig. |4j 

We measure Bk on ensembles at 4 different lattice spacings and a variety of pion and kaon 
masses. Renormalization is again performed nonperturbatively in the RI-MOM schemed For 
each lattice spacing, we interpolate the renormalized Bk to physical pion and kaon masses using 
various interpolators (see fig. [5]) and the resulting physical value is extrapolated to the continuum 
(see fig. [6])^] In addition to the finite volume corrections on pion masses, we also apply finite 
volume corrections to Bk^ 



a In fact, both the interpolation to the physical point and the continuum extrapolation are technically performed 
in one combined, global fit. 
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Figure 5: Interpolation of the renormalized lattice results for Bk to the physical pion and kaon masses. Note 
that the interpolation curves from different lattice spacings are almost on top of each other. 
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Figure 6: Continuum extrapolation of the renormalized lattice results for Bk at the physical pion and kaon 

masses. 



Both interpolation to physical pion and kaon masses as well as the continuum extrapolation 
turn out to be very mild. In addition, the effect of finite volume corrections is even smaller than 
it was on quark masses. Consequently, the systematic error of our result is less than half the 
statistical error and we obtain 

Bf- MOM (3.5GeV) = 0.5308(56)(23) (4) 

as our final, fully nonperturbative result. 

For further conversion of Q into other schemes, we use results for the 2-loop runningl 17 ! 18 ! 
Adding a conservative perturbative conversion uncertainty of 1%, we obtain 

B MS-NDR (2GeF) = o.5644(59) stat (25) sys (56) PT (5) 
B K = 0.7727(81) stat (34) sys (77) PT (6) 
The latter is compatible with the prediction Bk = 0.83 in 15 from a global CKM fit?^ 
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